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Abstract. - We prove the existence of low energy excitations in insulating systems at general 
filling factor under certain conditions, and discuss in which cases these may be identified as 
topological excitations. In the specific case of half-filling this proof provides a significantly 
shortened proof of the recent higher dimensional Lieb-Schultz-Mattis theorem. 



The classic 1961 result of Lieb, Schultz, and Mattis (LSM) [1], proving the existence of an 
excitation within energy ~ 1/L of the ground state for certain one dimensional spin chains, 
has had a large effect on the field. While it was then proven by Affleck and Lieb [2] that one 
dimensional systems either have gapless localized excitations or a local symmetry breaking, 
it has long been suspected that in higher dimensions there is a more interesting possibility of 
topological order [3]. 

One way to understand topological order is based on flux insertion. We give the physical 
argument here, and then discuss the difficulties in this argument which give rise to the need 
for the more careful argument of this paper. We consider a higher-dimensional system which 
is periodic in one direction. A spin-1/2 system can be mapped to a hard-core boson system on 
a lattice, with the presence or absence of a particle denoting spin up or down. If the particle 
system is superfluid, there is long range order in the x and y components of the spin in the 
original system, implying the existence of low energy excitations. On the other hand, if the 
particle system is insulating, it should be possible to insert 27r of gauge flux in the hopping 
of particles across a given line cutting the system, returning the Hamiltonian to the original 
one, but, for non- integer filling fraction, taking the system to an excited state which is very 
close in energy to the ground state. Using adiabatic flux insertion, this was suggested as a 
way to prove a higher dimensional LSM theorem [4]. The two possibilities would thus seem 
to be a superfluid system (or other system which resists flux insertion) which has low energy 
excitations, or an insulating system which has topological order. In either case, there is a 
state close in energy to the ground state. 

However, there is a serious problem with this argument. The definition of adiabatic flux 
insertion depends on the existence of a gap (but does not require anything about the magnitude 
of a gap); however, in a spin system with no disorder, any gap in the spectrum at zero flux 
must close at some non-zero value of the flux, thus making it impossible even to define an 
adiabatic flux insertion [5]. In a fractional quantum Hall system with disorder, there is a 
related problem. In a 1/3 quantum Hall state, the gap between the three approximately 
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degenerate ground states remains open for all values of the flux, even though the gap is 
exponentially small. However, this means that adiabatic flux insertion leaves the system in 
the lowest of the three states and does not produce the correct topologically excited state [6] . 
Thus flux must be inserted slowly, but not adiabatically. For fractional Hall systems, the 
experimentally relevant flux insertion is slow enough to avoid exciting states above the three 
almost degenerate ground states, but fast enough to "shoot through" the level crossing. 

A precise definition of this quasi-adiabatic flux insertion was given in [7] , and used to prove 
a version of the LSM theorem for spin systems with total S z — valid in arbitrary dimension. 
The flux was inserted in such a way that the ground state wave function was only disturbed 
near the flux insertion point. Physically, one can imagine that the flux was inserted sufficiently 
rapidly to prevent influences propagating around the system, but sufficiently slowly to avoid 
creating localized excitations. Compared to [7], the proof here is shorter and is generalized 
to systems at general filling fraction p instead of just spin systems with total S z = 0. It also 
explicitly constructs the low energy states for certain systems. Later, we discuss under what 
conditions these states can be identified as topological excitations. 

Definition of System — We consider systems defined on a lattice, using letters i,j, ... to 
denote lattice sites (throughout, the term "site" may also be used to refer to a unit cell 
comprised of several sites), and introducing a metric d(i,j) on the lattice. We assume that 
there is a conserved charge, Q = ^\ Qi, where Qi, the charge on site i, is quantized to be an 
integer (for example, Qi may be taken to be the number of particles on a site, or in a spin 
system of half-integer spins Qi may be equal to Sf + 1/2). We assume that the Hamiltonian, 
Ti can be written as H — ^i^- 1 , where the Ti 1 obey the finite range conditions [8,9]: (1) 
the commutator \H l ,0\ = for any operator O which acts only on sites j with d(i,j) > R 
where R is the range of the Hamiltonian; and (2) the operator norm ||7i 4 || < J for all i, for 
some constant J. Finally, the number of sites j with d(i,j) < R should be bounded by some 
number that we denote S. These finite range conditions include all short-range boson systems 
with a finite number of bosons allowed per site, as well as all short-range spin systems. 

We suppose that there are a total of V sites on the ci-dimensional lattice. We assume that 
the Hamiltonian is translationally invariant in one of the directions of the lattice, with period 
L: this is the length of the system. We define the filling fraction of the system to be p — Q/V. 
For each site i, we introduce coordinates (xi,yi) to specify sites, where Xi labels the distance 
along the length direction and yi labels the transverse directions. The coordinate x is periodic 
with period L and should be compatible with the metric d(i,j): any two sites i, j should have 
d(i,j) greater than or equal to the minimum over integers n of \xi — Xj — nL\. 

We also define a twisted Hamiltonian and a rotation operator. Let 

R(0)= J] expWQj - P)], (1) 

j,0<Xj<L/2 

where the product ranges over all sites j such that < Xj < L/2. We introduce two boundary 
condition twists as follows. Let Tig lt g 2 = X^^i e 2 ' wnere '^■g 1 e 2 ~ R(0i)W t R(—6i) if —R < 
Xl < R and U\ uH = R{-6 2 )WR{6 2 ) if L/2 - R < x t < L/2 + R and otherwise H\ u92 = W. 
Thus, angle 6\ defines a twist in boundary conditions between x = and x — 1, while angle 
62 defines a twist between x = L/2 and x = L/2 + 1. These twists correspond to inserting 
gauge flux between x = and x = 1 or between x = L/2 and x = L/2 + 1. Let there be a 
total of N c sites i with — R < Xi < R, and let WdgH^ 02 || < K for all i for some constant K. 
For a d dimensional system with linear size L and aspect ratio close to unity, N c ~ L d ^ 1 . 

Let ^0(^11^2) be the ground state of Hg lt g 2 , and let \t be the ground state of Ho,o = 
H. Note that U dlfi2 = R(-6 2 )Hg 1+g2fi R(9 2 ) and in particular Hg-e = R{6)HR{-6) and 
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$>o(9,— 8) = R(9)^!o. The second boundary condition angle thus appears to be redundant, 
since it is only the sum of the two angles, 8\ + 82, that affects the eigenspectrum of Ttg 1: g 2 . 
That is, the physical properties depend only on the total gauge flux 8\ + 82. However, the 
introduction of the second boundary angle will be extremely useful below. 

Results — The main result is that for a short-range Hamiltonian as defined above we are 
able to construct topologically excited states under certain assumptions on the spectrum of 
the Hamiltonian. The topologically excited states have a low energy but momentum that 
differs from the ground state momentum as given in Eq. (1 lloltill . We begin by specifying the 
conditions on the spectrum. 

Suppose that, for 81 = 8 2 = 0, there is a gap AE in energy between the ground state and 
a set of eigenstates, that we refer to as local excitations; we define Phigh to be the projection 
operator onto this space of states. There may be other eigenstates with energy below this 
gap; we define Pi ow = 1 — Phigh to project onto this space. We assume that for certain local 
operators O, namely O — dgHg.-g and O = dgR(9), the matrix elements of these operators, 
between ^0 and any other normalized state \P which is a linear combination of eigenstates 
below the gap and which is orthogonal to the ground state, are bounded by e| \0\ | for some small 
constant e. Physically, we are interested in this case because in many situations the matrix 
elements of local operators between the ground state and topologically excited states are very 
small, often exponentially small in the system size. Since the Hamiltonian is translationally 
invariant, the ground state, may be taken translationally invariant, with eigenvalue z : 
T^ = zq^q where T is the translation operator which increases the x-coordinate by 1. 

Then, we introduce the twist operator 

/>(p />oo 

Wi(<l>) = 6exp{- / d8 / dTeM-(rAE) 2 /{2q)}[w+ e { l T)-h.c.}}, (2) 
Jo Jo 

where the symbol denotes that the exponential is ^-ordered, in analogy to the usual time 
ordered or path ordered exponentials. We define Wi t g = dgHg y o and define e (ir) following 
[9]: for any operator A 

A(t) ee A(t)eM-(t^E) 2 /(2q)}; A ± (± l r) = ±- [ dt A(t)— 1— . (3) 

2n J ±lt + T 

The time evolution of operators is defined by A(t) = exp[ffl.i-gi\Aexp[— i7ii-$t], where we 
define the Hamiltonian Hi-.g — J2i o> w ith the sum ranging over —L/4 + R < Xi < L/4 — R. 
The Hermitian conjugate in Eq. @ of g {ir) is w^ g (— it), and Wtiff) is a unitary operator. 
We define W x = W^ir). 

To understand the definition for any operator A define A + (ir) to be the positive 
energy part of A at imaginary time it. That is, in a basis of eigenstates of the Hamilto- 
nian with energies JSf, we have Afj(ir) = Aij exp[— r(Ei — Ej)]Q(Ei — Ej), with the step 
function 0(x) = 1 for x > 0, 9(0) = 1/2, and Q(x) = for x < 0. The operator 
exp[— (tAE) 2 / (2q)]A + (ir) is a good approximation to A + (ir) in the following sense: the differ- 
ence \P high eM~(TAE) 2 /(2q)}A+{iT)y -P high A+(iT)y \is bounded by e X p[-(rAE) 2 /{2q)} exp[-q/2]\\A\\ 
for t < q/ AE and is bounded by exp[— tAE] \ \A\ \ for r > q/AE. This result follows from 
elementary integrations [9] and is a kind of uncertainty relation: A(t) is cutoff by the Gaussian 
at times of order 1/ AE and so can be used to approximate states with energies of order AE. 
In the limit q — > 00, when the approximation A + (ir) becomes exact, Eq. © adiabatically 
twists 8\ by 0; we instead keep q finite as chosen later. 

The use of r H.\fi to define the time evolution, rather than TLgfl as in [7] is a very useful 
technical trick which simplifies the proof. Using Hi-g, we have [Wi,0] = for any operator 
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O which acts only on a sites j with L/4 < < 3L/4. That is, W^i acts only on sites j with 
-L/4 < a;j < L/4. 

We define = Wf^o and define i?„ = |7Y| \I/ rl ) — (H) where the first and second pairs 
of angle brackets denote the inner product and ground state expectation value respectively. 
We prove below that 

E n < VJe n /2, (4) 

where e„ is given by Eq. (|14fl and 

K* w |r|¥„) - (T)exp[i2imp(V/L)]\ < e„. (5) 
For large L, we will find that e n is of order 

e„ ~ (K/AE)[N* exp{-cxLAE/S) + N^y/ Cl LAEexp(-L/^c) + N c ey/ Cl LAE], (6) 

where the constants Ci,£c are defined below. Thus, for a d-dimensional system with N c cx 
L d ^ 1 and with fixed AE and vanishingly small e, we have e n ,E n going to zero exponentially 
in L. 

One consequence of these results is a higher dimensional Lieb-Schultz-Mattis theorem for 
systems of with p(V/L) non-integer, a case which includes spin systems with half-integer spin 
per unit cell, total S z — 0, and odd width V/L [7] . If we set AE to be equal to the gap between 
the ground state and the first excited state, then e = as the ground state is the only state 
with energy less than AE. Then, if AE = klog(L)/c\L, the exponential exp(— ciLAE/8) 
decays as L _fc / 8 , and for large enough k the state "Ji will be less than AE in energy, while we 
can use the expectation value of T for state \&i to bound the overlap of Vfi with the ground 
state, leading to a contradiction and thus bounding AE. 

To show these results, we will need a few additional definitions defined in the next two 
sections. For use later, we introduce the operator 

W 2 (</>) = 0exp{- / d9 drexp[-(rAL) 2 /(2 9 )][w+ e (ir) - h.c.}}, (7) 
Jo Jo 

where u> 2 ,e = dgHo,e- We write Wi = Wz{2 , k). Here, the time evolution of operators is defined 
using the Hamiltonian Ti.2-e = J2i e> w ^ tn * ne sum ran ging over L/4 + i? < Xi < 3L/4 — R. 
Finally, we define W{4>) = 6exp{- d9 / °° dr exp[-(rAE) 2 /(2q)][wj(iT) - h.c.]}, where 
wg = dgli.0 t -g and where the time evolution is defined using the Hamiltonian TLg^-g. 

Twist in Boundary Condition — This and the next section show two basic facts about 
W,W\,W2- First, the difference |!4 7 (</>) , J'o — i?(0)\l/o| is exponentially small in q. Thus, 
W "twists" the boundary conditions, twisting 9\ by </> and 02 by ~<fi. Second, the difference 
1 1 W4 (0)^2(0) — W(^)|| is exponentially small in (c\LAE) 2 /q, and thus one can approximately 
factor the operator W into the product of two commuting operators Wi , W 2 ■ Physically, the 
operator W%(cp) twists 9% by while W2(<j>) twists 82 by —<\>. We will then be able to pick a 
q of order L such that both differences are exponentially small in L. Given these two results, 
Eqs. 1)10112(1 . it will be easy to derive the main results 1)415(1 above. The reader who would 
prefer to skip the derivation of these two results may skip to the section "Bound on Energy" . 

We start by using linear perturbation theory to compute dg^o(9, —9) = dgR(8)'i'o- We 
have Tig-gR(9)~$> = E R(9) , i' n . Taking derivatives of both sides with respect to 9 at 9 = 0, 
and working in a basis of eigenstates ^ of Ti. with energies JSj, we can compute the matrix 
element (d e R(9)) i0 if E t ? E Q : {d e R(9)) M = -{E % - E^- x {dgUg^g\ a . Thus, 

dg^ o (0, -9) = - V {d l He '~!} M ^ l + Z i0 *i = - r dT[{dgHg,-g) + {iT) - h.c.}V + Z* ,(8) 
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where we define the operator Z by Zij = (dgR(9))ij if £7, = Ej and Z^ = otherwise. The 
second equality in (JSJ) may be verified by computing the integral over r. In the limit q — > oo, 
the operator — J °° dT[(do'Tigi .-6>') + (* T ) — ^- c -] equals the integrand in the exponent of W and 
so in that limit, for a non-degenerate ground state, W(4>)^>o = R(4>)^a- 

We now consider the case of finite q and bound the difference |W / (0)^I'o — R(4>)^fo\. We 
use deTie-e = R(9)de>Tl$',-$'R(—9), taking the derivatives at 9' = 0, to rewrite W(</>) = 
6exp{- /* dO f™6Texp[-(rAE) 2 /(2q)]R(9)[w+(iT) - h.c.]R(-9)} , where the time evolu- 
tion of Wq is defined using the Hamiltonian H ,o = U. Thus, W(<p) = R((j>) exp{-0 / °° dr exp[-(rAE) 2 / (2q)][w+ {ir)- 
h.c] — (f>dgR(9)}, taking the derivatives at 9 = 0. 

Then, \W o - fl($)# | = \R(-<f>)W - * | and 

p OO 

ii(-0)W(^)*o = exp(-0 / dr{exp[-(rA£;) 2 /(2g)] ? I;+( l T)- W +^r)-/ l . C .}~0Z)*o. (9) 

Jo 

We now bound the norm of the state ^ = (— J °° dr{exp[— (rAi?) 2 /(2q , )]i2;Q"(ir) — u^J~(it) — 
h.c.} - Z)* - The norm |*| < |P Wgft *| + We have PhighZVo = 0, while using 

the bounds below Eq. © / °° dr{exp[-(rAi;) 2 /(2g)]w+(ir) - lu^(ir) - 7i.c.})* | < 

N c (K/AE)(2exp[-q] + exp[-g/2] v / 27n?). 

On the other hand, P iow ^ = Piow[~ r °° dr exp[~{rAE) 2 / {2q)]w^ {ir) - h.c. - d e R{9)}^ - 
The expectation values (dgR(9)) and (deHe^-e) both vanish, due to the translation symmetry 
of the ground state, so we may assume that Pi ow ^ is orthogonal to ^q. Then, we can 
use the bound on the matrix elements between and other states below the gap to show 
that \P low V\ < e^qN c {K/AE) + e| \d 9 >R(9')\ |. Therefore, |*| < N c (K/AE)[2exp(-q) + 
exp(-g/2) v ^] + ey/2^qN c (K/AE) + e\\d d 'R(9')\\. 

Therefore, \W(<f>)ty - R(<f>)^ \ < c 2 (4>), where we define 

c 2 {4>) = (b{N c (K/AE)[2exp(-q) (10) 
+ cxp(-g/2)V2^] + ey/2^N c (K/AE) + e\\d S 'R(9')\\}. 

Locality Bounds — The next result we need is a bound on ||Wi(</>)W2 (<?!>) — W(0)||. Note 
that W 1 ((j})W 2 ((j)) = 9exp{- f* 6.9 dTexp[-(rAE) 2 /(2q)][w£(iT) - h.c.]}, where w e = 
dgHg-g and where the time evolution is defined using the Hamiltonian 7ii ; #, +7^2;-e- This 
differs from the definition of W in only one respect: the use of Hi-fi, + Ti-2-,-9 to define the 
time evolution rather than Tlg^g. 

We next recall the finite group velocity result, proven in [7,8], that there exists a func- 
tion g(t,l), which depends on J, R, and the lattice structure, such that ||L4(t),_B(0)]|| < 
||^.||||-B||X}j3 , (*'(7')> wnere the sum ranges over sites j which appear in operator B, and 
lj = d(j, i) is the distance from j to the closest site % in the operator A. It was shown that 
there exists some constant c\ such for \t\ < c\l, g{c\l,l) is exponentially decaying in I for large 
I with correlation length £c- 

Now, we pick B = Tie-e — Hifi — 7~C2;-$, with < 2N C J. The operator B can be 

written as a sum of TC 1 , with each such Ti 1 having at most 5 sites. Then, we can bound the 
difference between the two different definitions of w~^{iT) using the two different definitions 
of time evolution. We have 

(27T)- 1 1 J 6t(it + T)- 1 exp[-(tAE) 2 /(2q)}{exp[iH8,-8t}w e exp[~iH8,-8t} - (11) 

exp[i(He,-e - B)t]w e exp[-i(H e -e - B)t]}\ 
<(2Tr)- 1 N c \\wg\\ f 6t(it + r)" 1 exp[-(tAE) 2 / (2q)} j dt'2 JSg{t', L/4 - 2R). 
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The integral fj dt'2JSg{t' , L/4-2R) is bounded by g{t,L/4) [7]. Thus, we find that the dif- 
ference using the two different definitions of time evolution is bounded by (2ir) ~ 1 N c 1 1 wg \ \ J dt(it- 
r) _1 exp[— (tAE) 2 /(2q)]g(t, L/4). To bound this integral, for t < ciL/4, we can use the lo- 
cality bound, while for t > c\L/4, we can use the Gaussian in t. The result gives a bound 
(2ir)- 1 2N c \\wg\\{g(c 1 L/4, L/4) + 2[ y ^/(AEc 1 L/4)}exp[~{c 1 LAE/4) 2 /(2q)}}. Integrating 
over r, and using \\wg\\ < N C K, we find | \W% ((f>)W2 (</>) — VF(</>)|| < c 3 (</>), where we define 



2-irq(K/AE)((t>/2w)2N 2 x (12) 




{g( Cl L/A, L/4) + 2 A ^'^ /4 exp[-(c 1 LA J B/4) 2 /(2g)]} 

Bound on Energy— We compute ^(Wf "WWf-W*). Note that Wx (-0) = wl{6). First 
we note that for sites i with L/4 < x % < 3L/4, (Wf "HWf - W) = 0, since W 1 commutes 
with H l for such sites. 

We now consider the change in expectation value of the energy for sites i with —L/4 < 
Xl < L/4. For such sites, \{W^ n WW?) = (W^W^WWfW^l Thus , we c an bound 
\{W^ n WW?) - (W)\ < 2J\W?W£y a - *o|- Then, we can combine Eqs. fTUTHl with the 
fact that i?(27rn)* = *o to bound \W^W^q - * | < c 2 + c 3 . Summing over all V/2 such 
i, we find 

£ | (^"HWf) - (W*) | = £„ < VJe n /2, (13) 

i 

where 

e n = 2[c B (27rn) + cs(27rn)]. (14) 
Finally, we pick q to be equal to Ci(L/2)AE to obtain the lowest possible E n , giving the 
exponential decay of e„ with L claimed in the results. 

We have thus bounded the difference in energies. The same argument would also work to 
bound the difference (Wi(— (j})H<pfiWi((j))} — (H) for general <j>. 

Expectation of Translation Operator — We next consider ( 1 I r „|T'|\I f ri ). We define a twisted 
translation operator, Tg$* = i?i(#i)i? 2 (#2)T, where = J^j Xj =i exp[z0i(Qj — p)] and 

^2(^2) = ]lj,xj=i/2+i exp[i0 2 (Qj - /?)]. Then, T^g, is a symmetry of H 9li e 2 and T e _ e = 
R(6)TR(-6), so that T0,_0\I>o(0, -0) = 20*0(6*, -#)• Since Q, is quantized, we have T 27m . = 
exp[-i27mp(V/Z,)]T. 

In this section, we will prove a bound on the difference \(Wi(— fyT^fiWi ((/))) Zq 1 — 1|. 
Taking = 27m, this will give a bound on the difference K^lT]^) — z exp[i2nnp(V/ L)}\ as 
desired. We have (W^-^T^qW^^Zq 1 = (^(-^^(^[TWi^T" 1 ]). Physically, this 
expectation value is close to unity for the following reason: the operator W\ inserts a twist 
between x = and x = 1, while TWi((f>)T~ 1 inserts the twist between x = 1 and a; = 2. The 
difference is the rotation of the sites with x = 1. Thus, Ri((j})[TWi((p)T^ 1 } is very close to 
the operator W\{<fi). We now make this argument precise. 

Note that (Wii-^Rxi^iTW^T- 1 }) = (W 2 (-(f>)W 1 (-(l))R 1 ((t>)[TW 1 ^)T- 1 }W 2 ((l))). 
Using previous results, we know that |*$W 2 (-0)Wi(-0) - #Ji?(-</>)| < c 2 (</>) + c 3 (</>). Now, 
consider [TWi(</!))T _1 ] W / 2 (0) v I'o. Let us define a new twisted Hamiltonian, 7i' e g 2 , in which 
the twist in boundary conditions is by angle 9\ between x = 1 and x = 2 and by angle 
2 between 1 = L/2 and a; = L/2 + 1. Specifically, H' 6u ^ = Rx{-6i)H' gug% Ri{Bx). De- 
fine W{(4>) = eexp{- d6 dr exp[-(TAE) 2 /(2q)][w'+ e (ir) - /i.e.]}, where ti/ 1)fl = cV^o 
and where the time evolution is defined using the Hamiltonian H! V g = THi-eT -1 . Cru- 
cially, = TW / i(0)T _1 . Similarly, define W in analogy to the definition of W, using 
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Then, following the same steps as above, we can show that |H^'(0)5 , o — i?i (— 0)i?(0)^E f o I < 
c 2 {4>) and \\W{(4>)W 2 (cP)-W'(<t>)\\ < c 3 (<£) so that {[TW^T-^Wi^o-R^-^R^ol < 
c 2 (0)+c 3 (0). Thus, \{W 2 {-<j,)W 1 {-<j 3 )R 1 {<j,)[TW 1 { ( t ) )T- 1 ]W 2 { ( t ) )) -1| < 2[c 2 (0) + c 3 (0)], and 
thus |(\I' n |T|\I' n ) — zo exp[i2Trnp(V/L)]\ < e n as desired. Note that in this derivation, the 
vanishing of the phase factors {d' g [Ri(—6')R(9')]) and Zoo = (d$'R(9')) is crucial; we chose 
the factor of Ilj,x J =i ex P[ — i@P] m the definition of R\ to make this phase factor vanish and it 
is this factor that led to the expectation value of T for state 

Discussion — The case of general filling fraction has been previously considered in one 
dimension, and used to study magnetization plateaus [10]. It was shown that gapped states 
could only exist at integer filling fraction, while at filling fraction p = p/q with p, q coprime 
there were found to be at least q low energy states. We have shown very similar behavior 
in higher dimensions for p(V/L) = p/q: if AE remains non- vanishing and e — > sufficiently 
rapidly as L — > oo then we find that the states \I/ n , for n = 0, q— 1 provide q degenerate low 
energy states. The nature of these q states is not known a priori. They may correspond to 
either discrete symmetry breaking or to topological order. If all local operators, not just the 
specific ones considered in the proof above, have matrix elements between the low-lying states 
which vanish as L — > oo, then we identify this as topological order. If some local operator 
has non-vanishing matrix elements between the local states, then there is long-range order in 
that operator [9], and we identify this as symmetry breaking. It is a limitation of the proof 
that the ability to construct low energy states depends on the width V/L of the system. We 
conjecture that if p = p/q then one can still find q low lying states for arbitrary, sufficiently 
large, V/L, but we are not able to prove this. 

The technique shown in this paper is very general, and can be used in any case in which one 
can identify a quantized conserved charge, including cases of spins transforming under higher 
SU(N) groups. Comparing to [11], the technique here creates a "global vison" excitation; if 
instead the flux is inserted along a finite line with endpoints, this technique can create a local 
"vison" excitation in two-dimensional spin systems as will be discussed elsewhere. 
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